A Measurement of the Branching Ratio of Kl — ► e + e 77 
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We report on a study of the decay Kl — > e + e~77 carried out as a part of the KTeV/E799 
experiment at Fermilab. The 1997 data yielded a sample of 1543 events, including an expected 
background of 56 ±8 events. An effective form factor was determined from the observed distribution 
of the e + e _ invariant mass. Using this form factor in the calculation of the detector acceptance, 
the branching ratio was measured to be B(Kl — > e + e~-fy, E* > 5MeV) = (5.84 ± 0.15 (stat) ± 
0.32 (sys)) x 10" 7 . 
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The process Kl — > e + e~77 occurs mainly through 
kaon Dalitz decay [Kl — > e + e~7) with an internally ra- 
diated photon. As such, measurements of the radia- 
tive decay can test radiative correction calculations and 
probe the form factor of the KlJJ* vertex. The study 
of Kl — > e + e~77 is also important because it is a signif- 
icant background to searches for Kl — ► 7r°e + e~, a mode 
which has an important direct CP- violating component 
and was a primary focus of the KTeV/E799 rare-decay 
program [jjj. 

The best existing measurements of the branch- 
ing ratio B(Kl — > e + e~j-f,E* > 5MeV) were made by 
the E799-I and NA31 experiments, which measured 
(6.5 ± 1.2 (stat) ± 0.6 (sys)) x 1CT 7 § and (8.0 ± 
1.5 (stat)^ 2 (sys)) x 10~ 7 ||, respectively, based on sig- 
nals of 58 and 40 events. This Letter presents a branch- 
ing ratio measurement with much higher precision, using 
data collected by the KTeV/E799 experiment at Fermi- 
lab in 1997. 

Figure [j] shows a plan view of the KTeV/E799 detec- 
tor. The Tevatron at Fermilab made an 800 GeV pro- 
ton beam, typically delivering 3.5 x 10 12 protons per 
minute. The beam was incident on a BeO target. Pho- 
tons from the target were converted by a lead absorber 
immediately downstream. Charged particles were then 
removed with magnetic sweeping. Collimators defined 
two neutral beams that entered the KTeV/E799 appa- 
ratus 94m downstream of the target. A 65 m vacuum 
decay region extended to the first drift chamber. The 
charged-particle spectrometer consisted of a dipole mag- 
net and four drift chambers (two upstream and two down- 
stream of the magnet) with ~ 100/xm position resolu- 
tion in both horizontal and vertical directions. The mag- 
netic field imparted a 205MeV/c horizontal momentum 
kick. The spectrometer had a momentum resolution of 
a{P)/P ~ 0.38% © 0.016%P, where P is in GeV/c. 

Photons were detected using an electromagnetic 
calorimeter, which consisted of 3100 pure Csl crystals, 
each 50 cm long 0] . Crystals in the central section of the 
calorimeter had a cross-sectional area of 2.5 x 2.5 cm 2 , 
and those in the outer region had a 5 x 5 cm 2 area. 
The calorimeter's energy resolution for photons was 
a(E)/E ~ 0.45% 2%/VW, where E is in GeV. The 
position resolution was about 1 mm. By comparing spec- 
trometer momentum to calorimeter energy, could be 
distinguished from . Nine photon veto assemblies 
(lead scintillator sandwiches) detected particles leaving 
the fiducial volume. 

Additional e ± separation was provided by eight 
transition radiation detector (TRD) stations located 
downstream of the last drift chambers J5j. Each TRD 
station consisted of a polypropylene felt radiator fol- 
lowed by two planes of multiwire proportional chamber 
(MWPC) filled with an 80/20% mixture of Xe/C0 2 gas. 
The pulse heights from the 16 MWPC planes, compared 



with the pion pulse height distributions from a sample 
of Kl — > it e*v decays, were used to calculate the confi- 
dence level that a given track was a pion. Pion rejection 
factors over 200:1 were possible with a 90% electron effi- 
ciency. 

The trigger system required hits in two scintillation 
hodoscopes just upstream of the Csl, as well as drift 
chamber hits consistent with two coincident charged par- 
ticles passing through the detector. The trigger used 
in this study counted the number of isolated clusters 
of in-time energy over 1 GeV with a special processor 
0; at least four such clusters were required. The to- 
tal energy deposited in the calorimeter had to be greater 
than 28 GeV. The scintillator hodoscope behind a lead 
wall downstream of the calorimeter vetoed events with 
hadronic showers in the calorimeter, and events with ac- 
tivity in the photon veto system were similarly discarded. 
Events which passed the hardware trigger requirements 
were read out into on-line CPUs, which performed ba- 
sic event reconstruction and applied a few loose event- 
topology and particle-identification cuts to select events 
to be written to tape. 

During offline event reconstruction, each event was 
required to contain two tracks from oppositely-charged 
particles originating from a common vertex. The recon- 
structed vertex had to be between 96 m and 158 m from 
the target. The tracks had to be at least 1 cm apart at the 
most upstream drift chamber, with an opening angle be- 
tween the tracks of at least 2.24 mrad. To select electrons, 
each track was required to point to a calorimeter cluster 
with an energy equal, within ±5%, to the track momen- 
tum measured by the spectrometer. The probability that 
a track was caused by a pion, formed by combining TRD 
chamber signals, was required to be less than 4% for each 
track. The total momentum of all decay products had to 
be under 216GeV/c and the total energy greater than 33 
GeV. Events with a total of four clusters were analyzed 
as Kl — * e + e~77 candidates. Events with a total of five 
clusters were analyzed as K — > tt q tt^ decays, where 77^, 
denotes a pion decaying to e + e~7. The latter mode was 
used both to check the simulation and to measure the 
total number of Kl decays in the data sample. 

In addition to the Kl e + e~jj signal mode, several 
potential background processes were considered. The 
Kl — > e + e~7 mode was a background when an addi- 
tional photon was present. This photon could have been 
radiated by an e while passing through matter in the 
detector ("external" radiation), or it could have been re- 
constructed from coincidental activity in the calorime- 
ter not associated with the decay ("accidental" activ- 
ity). K — > tt^tt^j was a background when one photon 
was missed. Similarly, Kl — > 7r°7r 7T£, could have been 
a background but was suppressed because missing three 
photons generally resulted in a reconstructed mass well 
below the Kl mass. Background from Kl — > ir^e^v plus 
two additional accidental photons with the 7r ± misiden- 
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tified as e by the calorimeter, was eliminated using the 
TRDs. 

Monte Carlo (MC) simulations of the detector 
were performed to calculate the acceptances for 
Kl — > e + e _ 77 and K — > K Q Tt D and the misidentification 
rates for Kl — > e + e~7, if — > tt ^ and AT^ — * 7r ± e =F ^. 
The effect of accidental clusters was simulated by 
overlaying on the MC events with data events taken 
with a random trigger that had a rate proportional 
to the beam intensity. The finite angle of external 
radiation was simulated. Simulations of Dalitz de- 
cays included radiative corrections to Oia 2 EM ) based 
on the work of Mikaelian and Smith 0. The ra- 
diative correction simulation also quantitatively pre- 
dicted the ratios T(Kl — > e + e~-f)/T(KL — > 77) and 
T(Kl — > e + e _ 77)/r(A'£ — > 77) and kinematic distribu- 
tions such as photon energy in the Kl center-of-mass, 
E*. The simulations of Kl — > e + e~"f"f and Kl — * e + e~7 
used the Bergstrom, Masso, and Singer (BMS) form fac- 
tor which includes a parameter ax* describing the 
relative strengths of the intermediate vector and pseu- 
doscalar meson amplitudes in the Kl — > 77* vertex. 

The number of Kl decays was measured by the 
K — > 7r°7T£, decay mode, as in the Kl — > ir°e + e~ anal- 
ysis 0. There were (2.642 ±0.015 (stat) ± 0.024 (sys) ± 
0.091 (BR)) x 10 11 decays in the fiducial region, where 
the third term is from uncertainties on the branching ra- 
tios of K L -> vr°7r° and tt° -> e+e-7 @. 

Several cuts on reconstructed quantities were made on 
events with exactly four calorimeter clusters to identify 
Kl — > e + e~77 signal candidates. First, to ensure that 
candidate events did not include radiated photons below 
the nominal infrared cutoff of 5 MeV, E* was calculated 
for each photon and was required to be at least 8 MeV (to 
allow for finite detector resolution) . Since the MC sample 
was subjected to the same requirement, the acceptance 
correction procedure yielded a branching ratio valid for 
E* > 5 MeV, allowing direct comparison to theoretical 
predictions as well as to other published measurements. 

In order to reduce backgrounds in which one or more 
particles are missing or mismeasurcd in the detector, or 
which involve accidental activity, a momentum balance 
cut was imposed. The square of the component of the to- 
tal momentum of the decay perpendicular to a line drawn 
from the target to the vertex was required to be less than 
300 (MeV/c) 2 for K L -> e+e~77 candidates. 

The invariant mass of the decay, m ee77 , was calcu- 
lated assuming that the tracks were electrons. This mass 
was required to be within 11 MeV/c 2 of the neutral- kaon 
mass. Figure |^ shows the distributions of m ee77 , after 
all other cuts, for data, signal, and the two main back- 
grounds. Values of m ee77 for K — > vr ^ decays tended 
to be less than the kaon mass because one photon was 
missing. The excess in data at very low m ee77 was from 
Kl — > 7r^7r^7r^j decays with three missing photons. The 
Kl — > e+e _ 7 decays often had a mass above the kaon 



mass because most such decays had an extra, accidental 
photon. The peak at the kaon mass in the Kl — ► e + e~7 
distribution is caused by events with external radiation. 

The distributions of minimum angle, 6 m i n , between 
any photon and any electron in the center of mass of the 
decay are plotted in Fig. || after all other cuts. The small 
angle events were dominated by Kl — * e + e~j with an 
externally radiated photon. The large angle events were 
dominated by K — > 7r°7r^ , where the electron directions 
were relatively uncorrelated with photon directions. Sig- 
nal events were required to have 0.05 < 9 m i n < 0.63 rad. 

There were 1,543 Kl — > e+e~77 candidate events that 
satisfied all requirements. The distribution of e + e~ mass 
in these events after background subtraction was used to 
determine the form factor. However, since the MC did 
not include radiative corrections to Kl — * e + e~77, an 
effective parameter, ckes, was measured rather than the 
true a if* . We simulated the effect of different effective 
form factors in the MC by reweighting Kl — * e + e~77 
and Kl — * e + e _ 7 MC events with a ratio of the BMS 
form factor squared over the generated form factor 
(where a^s = —0.28) squared. The best fit to the data 
was found at cues — +0.016 with x 2 = 4.3 for 9 degrees 
of freedom. Fig. |] shows the ratio of data to MC for this 
value and for a^s = —0.30. 

The statistical uncertainty was found by using 22 sub- 
sets of MC events, each the same size as the data. The 
spread of different cues found by using each subset in 
place of the data was taken as the statistical uncertainty, 
0.083. A sample of 12 million Kl — * Tr°TT TT° D decays was 
used to check the MC simulation of the acceptance as a 
function of m ee . A ratio of the data over MC results in a 
slope of 0.038±0.096 (GeV/c 2 )" 1 and a ±la range of this 
slope was fed back into the Kl — * e + e~77 simulation, re- 
sulting in a systematic uncertainty of ±0.032 on aEff. 
Uncertainty in our ability to reconstruct events with 
Kl — > 7r 7r 7r^, decays, where the two (usually closely 
spaced) tracks both passed through one of the beam re- 
gions of the first drift chamber led us to assign an uncer- 
tainty to aEff of ±0.027. Other systematic uncertainties, 
similar to those assessed below for the branching ratio, 
were evaluated, but the results were negligible. These 
systematic uncertainties are combined to give a total un- 
certainty of ±0.042 on aEff. 

The acceptance for Kl — > e + e~77 with generated 
E* > 5 MeV was found to be 0.984% using MC 
reweighted to have a^s — +0.016. A Monte Carlo 
simulation that was normalized to the number of Kl 
decays in the detector predicted 31 background events 
from Kl — > e + e~j (and Kl — > e + e~77 with true E* < 
5 MeV) and 25 background events from K —> n ^- 
Thus, the branching ratio of K L —> e + e~77, E* > 5 MeV 
is (5.84 + 0.15 (stat)) x 10~ 7 . 

Several sources of systematic uncertainty for the 
branching ratio measurement were considered. The nor- 
malization calculation contributed uncertainties as men- 
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tioned previously. A systematic uncertainty of ±0.49% 
was assigned in order for the MC to match the relative 
excess of simulated events at small angles (Fig. ||) . An- 
other uncertainty was assigned to account for uncertainty 
in the inefficiency for tracks in the neutral-beam region of 
the drift chambers (the beam region was 0.965 ±0.028 ef- 
ficient), by reweighting MC events (signal, backgrounds, 
and normalization) based on the number of tracks in the 
beam region in each event; the resulting 0.49% drop in 
branching ratio was taken as a symmetric uncertainty. 
Varying aEft by our combined uncertainty of ±0.093 
changed the acceptance and Kl — * e + e~j background, 
and the resulting ±2.25% change in branching ratio was 
taken as an uncertainty. Finally, uncertainties due to 
inaccuracies in the MC simulation of the detector were 
estimated by varying one cut at a time over a reasonable 
range simultaneously for data, signal, backgrounds, and 
normalization. The most significant of these uncertain- 
ties were from the transverse- momentum cut, the cut on 
the opening angles between the the two tracks in the lab 
frame, and the m ee77 cut. When added in quadrature, 
all of the cut variations contributed a systematic uncer- 
tainty of ±3.92%. The combined systematic uncertainty 
on B{K L -> e+e~jj,E* > 5MeV) was ±5.44%. 

MC numerical integration of the tree-level QED partial 
width, using the BMS form factor with a^* = +0.016 
and B{K L -> 77) = (5.86 ± 0.15) x 10~ 4 |], predicts 
B(K L 



' + e~jj,E* >5MeV) = (5.70 ± 0.14) x 10" 



based on B(Kl 



e + e 7) 



(9.1 ±0.5) x 10~ b , in agree- 



ment with the measurement. 

In summary, we have determined the branching ratio 
of Kl — > e + e~77, with an infrared threshold of E* > 
5MeV, to be (5.84 ± 0.15 (stat) ± 0.32 (sys)) x 10~ 7 . 
This calculation uses an effective form factor, which 
does not include effects from radiative corrections, with 
a Eff = +0.016 ± 0.083 (stat) ± 0.042 (sys). This mea- 
surement is an improvement in precision over previous 
measurements, and agrees with the predicted value at 
tree level. 
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FIG. 1. The KTeV/E799 detector configuration for rare 
decay studies. 
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FIG. 2. The e + e~77 invariant mass distributions for the 
data, the Kl — > e + e~ r y r y MC, and the background MC sam- 
ples. The MC distributions use the normalization calculation 
and predicted branching ratios (for Kl — > e + e _ 7 ( 7 )) or PDG 
branching ratio (for K — > tt°iy%). 
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FIG. 4. The raw e + e~ invariant mass distribution for the 
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ant mass(m ee ) distributions, for MC using «BfF = 0.016 and 
ttEft = —0.30. The error bars are from data statistics. The 
MC distributions were normalized to have the same integral 
as the data distribution. 
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